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Abstract 

We study the implications of the global U(1)r symmetry present in minimal 
lepton flavor violating implementations of the seesaw mechanism for neutrino masses. 
In the context of minimal type I seesaw scenarios with a slightly broken U(1)r, we 
show that, depending on the i?-charge assignments, two classes of generic models can 
be identified. Models where the right-handed neutrino masses and the lepton number 
breaking scale are decoupled, and models where the parameters that slightly break 
the U(1)r induce a suppression in the light neutrino mass matrix. We show that 
within the first class of models, contributions of right-handed neutrinos to charged 
lepton flavor violating processes are severely suppressed. Within the second class 
of models we study the charged lepton flavor violating phenomenology in detail, 
focusing on /i — > &y, /i — > 3e and /z — e conversion in nuclei. We show that sizable 
contributions to these processes are naturally obtained for right-handed neutrino 
masses at the TeV scale. We then discuss the interplay with the effects of the right- 
handed neutrino interactions on primordial B — L asymmetries, finding that sizable 
right-handed neutrino contributions to charged lepton flavor violating processes are 
incompatible with the requirement of generating (or even preserving preexisting) 
B — L asymmetries consistent with the observed baryon asymmetry of the Universe. 
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1 Introduction 



The observation of neutrino flavor oscillations constitutes an experimental proof of 
lepton flavor violation [1]. In principle, other manifestations of such effects could 
be expected to show up in the charged lepton sector as well. However, the lack of 
a definitive model for neutrino mass generation implies that conclusive predictions 
for lepton flavor violating processes can not be made, and even assuming a concrete 
model realization for neutrino masses, predictions for such effects can only be done 
if the flavor structure of the corresponding realization is specified. 

A problem that one faces when dealing with charged lepton flavor violating phe- 
nomenology is related with the arbitrariness of the free parameters that define these 
observables. In this regards the minimal lepton flavor violation (MLFV) hypothesis 
[5, 6, 7, 8] is a very useful guide for constructing predictive models in which lepton 
violating signals are entirely determined by low-energy neutrino data. The relation- 
ship between the free parameters and neutrino observables can arise via either a 
restrictive MLFV hypothesis or in some cases due to the intrinsic structure of the 
corresponding model (whenever the number of free parameters is comparable to the 
number of neutrino observables). 

In the canonical seesaw (type-I) the kinetic sector of the model is invariant under 
a global large flavor symmetry group Gp = SU(3) e x SU(3)i x SU(3)n, where £, 
e and N denote triplets in flavor space constructed from the electroweak lepton 
doublets and singlets and RH neutrinos. In addition there is an invariance under 
an extra global U(1)r [4, 8]. In principle the charges associated with this global 
transformation (hereafter denoted by R) are arbitrary, and thus different .R-charge 
assignments define different models with their own consequences for lepton flavor 
violating phenomenology. In particular, models for which the R charges allow for 
large Yukawa couplings and TeV RH neutrino masses should lead to sizable charged 
lepton flavor violating processes. 

Consistent models of O (TeV) RH states and large Yukawa couplings are achiev- 
able if cancellations among different pieces of the light neutrino mass matrix are 
allowed, and the RH neutrino mass spectrum is not strongly hierarchical [9] 1 . The 
class of TeV scale seesaw models arising from the presence of the U(1)r symmetry 
are expected to be in that sense different: no cancellations are needed because the 
suppression in the effective light neutrino mass matrix is no longer constrained to 
be related with the heaviness of the RH neutrinos. 

In this paper we study the implications of a slightly broken U(1)r symmetry in 
minimal type-I seesaw models (with two RH neutrinos) 2 3 . We show that the intrin- 
sic structure of the relevant models leads to a flavor pattern completely determined 

■"■In the case of a large hierarchy among the different RH neutrino masses the one- loop finite 
corrections to the light neutrino mass matrix can exceed the corresponding tree-level contributions. 
Neglecting such corrections can in this case lead to a model inconsistent with neutrino data [10]. 

2 Models with an arbitrary number of RH neutrinos were considered for the first time in [11]. 

3 For related discussion within Type III and mixed Type I+III seesaw scenarios c.f. [12]. 
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by low-energy neutrino observables, thus realizing in that way the MLFV hypoth- 
esis. This feature in addition to a slightly broken U(1)r leads to sizable \i — > ej, 
fi — > 3e and // — e processes. This is in contrast to models where U(1)l is slightly 
broken so that lepton number and lepton flavor violation occur at the same scale. 

In models with slightly broken lepton number it has been shown that leptogenesis 
is reconcilable with large charged lepton flavor violating rates, as the washouts 
induced by the RH neutrino states are controlled by the amount of lepton number 
violation [13]. We show however, that in the class of models considered here this 
statement does not hold, and that indeed leptogenesis and sizable charged lepton 
flavor violating effects are mutually exclusive. However, since both phenomena cover 
non-overlapped regions of parameter space their analyses are complementary. We 
therefore also explore the constraints on these models derived from the requirement 
of not erasing — via the RH neutrino dynamics — a preexisting B — L asymmetry 
below the value consistent with the observed baryon asymmetry of the Universe. 

The rest of this paper is organized as follows: in section 2 we discuss in detail the 
scenarios arising from the different -R-charge assignments. In section 3 we discuss 
the phenomenology of charged lepton flavor violating decays in the representative 
models. In section 4 we analyze the implications of such constructions for scenarios 
of high scale baryogenesis by quantifying — via the Boltzmann equations — the effects 
of the RH neutrino dynamics on preexisting B — L asymmetries. In section 5 we 
present our conclusions and final remarks. Explicit formulas used in the calculation 
of the different lepton flavor violating processes under study are given in appendix A. 

2 The setups 

The kinetic and gauge interaction Lagrangian of the standard model extended with 
two RH neutrinos exhibits a global G = U(3) e X U{3)txU{2)^ symmetry. Factorizing 
three U(l) factors from G, the global symmetry can be rewritten as{7(l)yx{7(l)z,x 
Z7 (l)ij x Gp where U(1)y,l can be identified with global hypercharge and lepton 
number whereas the U(1)r is a "new" global symmetry [4, 8]. The remaining direct 
product group Gp = SU(3) e x SU(3)e x SU(2)n determines the flavor symmetry 
which is explicitly broken in the Yukawa sector. 

In minimal lepton flavor violating seesaw models the Yukawa (mass) matrices are 
treated as spurion fields transforming under Gp in such a way that the corresponding 
Yukawa (mass) terms in the leptonic Lagrangian remain invariant under the global 
flavor symmetry. The usual procedure is then based on an effective theory approach 
in which a set of non-renormalizable effective operators are constructed from the 
spurions. With the operators at hand and under certain restrictions on Gf the 
lepton flavor violating effects can be estimated by means of low-energy neutrino 
data [5, 8]. Here, instead, we explicitly consider the seesaw Lagrangian with a 
slightly broken U(1)r and classify the possible realizations according to the -R-charge 

4 An exception are the explicit MLFV models discussed in Ref. [7]. 
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assignments. Under this consideration, in the models featuring sizable lepton flavor 
violating effects, the flavor structure is determined by low-energy observables as well 
(up to a global normalization factor), not as a consequence of a restricted MLFV 
hypothesis but by the intrinsic structure of the resulting models. 

Depending on the .R-charge assignments two classes of generic models can be 
identified. Let us discuss this in more detail. Requiring U(1)r invariance of the 
charged lepton Yukawa terms determines R(e) in terms of R(£,H). After fixing 
R{H) = 0, to avoid charging the quark sector, the remaining charges can be fixed 
by starting with i?(iV lj2 ). In order to have sizable lepton flavor violating effects 
both the iV lj2 Majorana mass terms should be suppressed by i?-breaking parameters 
(generically denoted by e), so R(Ni t2 ) 7^ 0. Thus one has only three possibilities: 



(A) R(Ni) = R(N 2 ), (B) R(N ± ) = -R(N 2 ) and (C) |i2(iVi)| ^ \R(N 2 )\. The 



phenomenology of case (C), however, is expected to be similar to the one arising 
from models with .R-charge assignments of type (A). The reason being that in that 
case the Ni—N 2 mixing is always U(1)r suppressed, which in turn implies suppressed 
charged lepton flavor violating processes (equivalently, the effective neutrino mass 
matrix will not be e suppressed, forcing tiny Yukawas or heavy N\ i2 states). This 
is to be compared with models based on .R-charge assignments of type (B), where 
the Ni — N 2 turns out to be maximal and a set of unsuppressed £ — N Yukawa 
couplings can be always obtained by properly chosing R(£) (in these models the 
effective neutrino mass matrix involves extra e suppression factors). 

In summary depending on the -R-charge assignments two classes of generic models 
can be identified: models in which the mechanism that suppresses the light neutrino 
masses propagates to lepton flavor violating observables, thus rendering their values 
far below planned experimental sensitivities; and models in which the mechanism 
decouples in such a way that lepton flavor violating effects become sizable. In what 
follows, in order to illustrate this is actually the case, we will discuss two examples 
of models of type A and B. 

Type A models [8]: R(N a ) = +1 and R(ii,e u H) = (H being the Higgs elec- 
troweak doublet) 

In this case, in the basis in which the charged lepton Yukawa couplings and the 
Majorana RH neutrino mass matrix are diagonal, the Lagrangian reads 



diag(Y/Vi, Yn 2 ) are the Yukawa coupling matrices (we denote matrices in bold-face). 
The dimensionless parameter e 1 slightly breaks U(1)r whereas, due to the 
assignment L(N) = 1, the lepton number U(1)l factor is broken by fi. With this 
setup the 5x5 neutral fermion mass matrix can be written as 




(1) 



Here H 



icr 2 H*, C is the charge conjugation operator and Y e , X* and Vjv 




(2) 
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where (H) = v ~ 174 GeV. In the seesaw limit, which in this case reads vX <C efiY^ 
the effective light neutrino mass matrix is given by 

r 0=1,2 iVa 

with X a = (Aea> A^a, A ra ). The corresponding light neutrino masses are obtained 
from the leptonic mixing matrix U = V P (with V having a CKM form and P = 
diag(e J9i , e 1 ^, 1) containing the Majorana CP phase) after diagonalization: 

U T mfU = mf. (4) 

In the 2 RH neutrino mass model the constrained parameter space enforces one of 
the light neutrinos to be massless. Thus, in the normal hierarchical mass spectrum 
case m vi = and m U2 < m U3 whereas in the inverted case m U3 = and m Vl < m U2 . 

Since the dimension five effective operator is U(1)r invariant the neutrino mass 
matrix does not depend on e; the suppression ensuring light neutrino masses is solely 
provided by the lepton number breaking parameter /1. On the other hand, the RH 
neutrino mass spectrum is determined by 

M N = e 2 /iYjv . (5) 

From this expression it can be seen that as long as the U(1)r global symmetry is 
an approximate symmetry of the Lagrangian (e< 1) the RH neutrino mass scale is 
decoupled from the lepton number violating scale. Thus, the RH neutrino masses 
do not lie at the same scale at which lepton number breaking takes place. 

Assuming Yn, X < 0(1) an estimation of the lepton number breaking parameter 
jj, ~ 10 15 GeV can be obtained using a/ Am^ ~ 0.05 eV [14] as a measure of the 
largest light neutrino mass in these scenarios. From this estimation and eq. (5) it 
can be seen that values of e of the order of ~ 10~ 6 allow to lower the lightest RH 
neutrino mass below 1 TeV. 

Formal invariance of the Lagrangian under Gf is guaranteed if the Yukawa ma- 
trices, promoted to spurion fields, transform according to 

Y e ~ (3 e ,3^, 1 N ) , A* ~ (l e ,3£,2jv) , Y N ~ (l e , 1^,3^) . (6) 

The constraints imposed by Gp imply 

A = X e ® X N , (7) 

where A^ is a SU(3)g triplet and Ajv a SU(2) N doublet in flavor space. Accordingly, 
in these kind of models a unequivocal determination of the flavor structure via the 
MLFV hypothesis is possible by means of a restrictive flavor symmetry G' F C Gf- 
Though several possibilities may be envisaged we do not discuss further details since, 
as we show below, in this type of models contributions to lepton flavor violating 
processes of charged leptons are always negligible. 
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Type B models : R{N 1 ,l h e i ) = +1, R(N 2 ) = -1, = 0. Changing R charges 

to L charges, this case resembles models where lepton number is slightly broken (see 
for example [7, 15, 16, 17]). The Lagrangian is given by 

£ = lY e eH - / Af N X H - e x I A 2 * A,// - |jVf CM N 2 - y.y.Yj CM aa iV a + h.c. . 

(8) 

The eA,iv parameters determine the amount of U(1)r breaking and are thus tiny. 
The diagonalization of the Majorana RH neutrino mass matrix leads to two quasi- 
degenerate states with masses given by 

M n + M 22 

M Nl 2 = M =f e^v . (9) 

In the basis in which the RH neutrino mass matrix is diagonal the Yukawa couplings 
become 

A fc a -> [Afei + (-l) a e A A fe2 ] , (k — e, /i,t and a = 1, 2) , (10) 

and the 5x5 neutral fermion mass matrix is similar as in type A models. However, 
due to the structure of the Yukawa couplings the effective light neutrino matrix, up 
to C(ejye^), has the following form 

mf =-^|Ai||A|(Ai*® A* +A*®v) , (11) 

with 

Mii + Moo ex - 

A =A2 AM~ e~ N Xl ■ (12) 

Here with the purpose of relating the flavor structure of these models with low 
energy observables, and following ref. [7], we expressed the parameter space vectors 
Ai, A in the light neutrino mass matrix in terms of their moduli |Ai|, |A| and unitary 
vectors Ai,A. 

Note that in these models lepton number is broken even when U{1)r is an exact 
symmetry of the Lagrangian. However due to the Yukawa structure and degeneracy 
of the RH neutrino mass spectrum at this stage m® ff = 0. Although a non-zero 
Majorana neutrino mass matrix arises only once the R breaking terms are present 
this does not imply that in the absence of lepton number violating interactions a 
Majorana mass matrix can be built. In that case — as can seen from eq. (8) — only 
Dirac masses can be generated, as it must be. 

Since e\ 1 small neutrino masses do not require heavy RH neutrinos or small 
Yukawa couplings, thus potentially implying large lepton flavor violating effects. In 
that sense, as already stressed, these models resemble those in which lepton number 
is slightly broken but with lepton number as well as lepton flavor violation taking 
place at the same scale M. 
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In contrast to models of type A, in this case due to the structure of the light 
Majorana neutrino mass matrix the vectors Ai and A can be entirely determined 
by means of the solar and atmospheric mass scales and mixing angles, up to the 
factors |Ai| and |A|, without invoking a restrictive MLFV hypothesis. The relations 
are different for normal and inverted light neutrino mass spectra [7]: 

• Normal hierarchical mass spectrum 

Aj = |Aii Ai = ^ (yr+pw + v / r 3 p£/2*) , (13) 

A = |A| A = ^ (y/lVpW - v / l 3 P#2*) , (14) 
where Ui denote the columns of the leptonic mixing matrix and 



a/1 + r - a/t 
a/1 + r + y/r 

Inverted hierarchical mass spectrum 



2 

m. 



P = V ^— V F , r= 2 "^ 2 - . (15) 
a/1 + r + y/r m u 3 - m t 2 



Ax = |Ai| A x = ^ {s/l+p~W + a/W^i*) , (16) 
A = |A| A = ^ [y/l + ~ P W - TWtfij , (17) 



with 



y/l+r - 1 mL - m 2 



P= nF= — » r = • (1 



\/l + r + 1 ' m 



With these results at hand we are now in a position to calculate the most relevant 
lepton flavor violating processes, which we discuss in turn. 



3 Lepton flavor violating processes 

In type A models the RH neutrino masses can be readily at the TeV scale for 
e ~ 1CT 6 . Since the Yukawa couplings scale with e as well, type A models are — in 
that sense — on the same footing as the canonical type-I seesaw model i.e. TeV RH 
neutrino masses imply tiny Yukawa couplings and thus negligible charged lepton 
flavor violating decay branching ratios. The main difference is that in the canonical 
case, sizable charged lepton flavor violation can still be induced via fine-tuned can- 
celations in the effective neutrino mass matrix, while no such effects are possible in 
the minimal type A models, since the relevant couplings are completely determined 
by light neutrino mass and mixing parameters. 
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Type B models, in contrast, may exhibit naturally large Yukawa couplings even 
for TeV scale RH neutrino masses (or even lighter masses, depending on the value 
of the U(1)r breaking parameter e\). Accordingly, several charged lepton flavor 
violating transition rates — induced by the RH neutrinos at the 1-loop level — can 
in principle reach observable levels. In what follows we study the allowed mass 
and Yukawa normalization factor ranges by considering the following lepton flavor 
violating processes: k — > lj-y, U — > 3lj and [/, — e conversion in nuclei. 



3.1 li — > Ijj processes 

Among these lepton flavor violating processes, presently the \i — > e^y transition is 
most severely constrained. The MEG collaboration recently established an upper 
bound of 2.4 x 10 -12 at the 90% C.L. [18]. For r — > and r — > /ry on the other 
hand, the bounds are 3.3 x 1(T 8 and 4.4 x 1(T 8 at 90% C.L. [20], respectively. In 
the limit mi- <C m;. the partial decay width for Zj — > processes reads [21] 



256 7T 2 M w 1 7 1 1024 7T 4 Myy 



where Myy is the W ± mass, aw = 5 ,2 /4vr and the elements of the diagonal matrix 
G 7 are given in eq. (37) in the appendix. This function is such that in the limit 
M Na 3> M w , (G 7 ) aa — > .\ /{"j - / 2.\ / y . The corresponding decay branching ratios are 
determined from the partial decay width after normalizing to T^ ot = Ht^, with n t 
the U charged lepton mean lifetime. In the limit r a 1, using eq. (5) and taking 
into account the Yukawa rescaling A — > e\ the decay branching ratio in type A 
models can be written as 



5 

BR(/j -> Li) ~ — %—r ^% -j- (X Y N - 2 \n ..' . (20) 
Kl JU 4096vr V 4 e 4 r^ ot V JlJ 



Thus, assuming 0(\,Yff) ~ 1, for which \i ~ 10 15 GeV and taking e = 10~ 6 , the 
value required for C(TeV) RH neutrino masses, we get BR(/i — > cy) ~ 10 -30 . This 
behavior, being extensible to other lepton flavor violating processes, shows that in 
type A models lepton flavor violating effects are negligibly small. 

In type B models in contrast such lepton flavor violating effects may be sizable. 
Using expression (10) for the Yukawa couplings, neglecting the piece proportional 
to e\ and taking the limit M Na 3> M w the decay branching ratios can be expressed 
in terms of the parameters Ai: 



a m 5 |Ai 



2 



BR(/j la) ~ — - -rjj Xn X* . (21) 

v 3U 1024vr 4 M 4 rl^ 3 

Since the components of the unitary vector Ai are entirely determined by low-energy 
observables (see eqs. (13) and (16)) the size of these branching ratios — and all the 
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M N [GeW] M w [GeV] 

Figure 1: Decay branching ratio BR(n — > cy) normalized to |Ai| 4 for normal 
light neutrino mass spectrum as a function of the common RH neutrino mass (left 
hand side plot). The upper horizontal dashed line indicates the current limit on 
BR(fi ej) from the MEG experiment [18], whereas the lower dotted one marks 
prospective future experimental sensitivities [19]. The corresponding bounds on |Ai| 
from the present [18] (upper gray band) and prospective future [19] (lower brown 
band) experimental searches are shown in the right hand side plot. The widths of 
the bands are due to the uncertainties in the neutrino mass matrix parameters. The 
results for the inverted neutrino spectrum are very similar and are thus not shown 
separately. 

others discussed below — are controlled only by the parameters M and |Ai|, thus 
implying that for sufficiently light RH neutrino masses and large |Ai| these processes 
may be measurable. 

In order to quantify the size of these lepton flavor violating effects we ran- 
domly generate neutrino masses, mixing angles and Dirac and Majorana CP vi- 
olating phases in their 2a ranges for both normal and inverted hierarchical light 
neutrino mass spectra [1]. We also randomly generate the parameters |Ai| and M in 
the ranges [10~ 5 , 1] and [10 2 , 10 6 ] GeV allowing RH neutrino mass splittings in the 
range [1CT 8 , 10~ 6 ] GeV. With the numerical output we calculate the different lepton 
flavor violating decay branching ratios from eq. (19), using the full loop function 
given in the appendix, eq. (37). 

We find that radiative r decay rates are always below their current bounds and 
barely reach values of 10~ 9 for RH neutrino masses around 100-200 GeV (values 
exceeding the current bounds are not consistent with the seesaw condition, that for 
concreteness we take as moM/v 1 < 10 _1 ), we thus focus on the \i — > process. 
The results for the normal mass spectrum case are displayed in fig. 1 as a function 
of the common RH neutrino mass, = M. We observe that BR(/i — > cy) can 
reach the current experimental limit reported by the MEG experiment [18] for RH 
neutrino masses M N < 0.1 TeV, 1 TeV, 10 TeV provided |Ai| > 2 x 10^ 2 , 10~\ 1, 
respectively. The results for the inverted light neutrino mass hierarchy are very 
similar and consistent with these values. Finally we note that the widths of the 
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bands in fig. 1 (and similarly for all the other considered processes below) are solely 
due to the uncertainties in the light neutrino mass matrix parameters (mainly #13 
and the CP violating phases) and can thus be improved with more precise light 
neutrino data. 



3.2 17 



processes 



The decay branching ratios for these processes have been calculated in [21] . The most 
constrained process in this case is \x~ — > e + e~e~ for which the SINDRUM experiment 
has placed a bound on the decay branching ratio of 10 -12 at 90% CL. [22]. For 
t~ e + e~e~ and t~ — > fi + fi~fi~ the current bounds are 2.7 x 10~ 8 and 2.1 x 
10 -8 , respectively [23]. The decay branching ratios for these lepton flavor violating 
reactions are given by [21] 

mi, ' ' 



a 



w 



24576tt 3 ML rh 



W 1 Total 



-^Box + ** 



+4s 



rplilj 
r Z 



2s 



rplilj 

c z 



TpUlj 

7 



+ 16s w Re 



piiij 

7 



Z 



Box 



-48s^Re 



+ 32s i w \G l ^\ 2 



G l * lj * 
In— I- 



11 

T 



(221 



Here sw = sinfljy, where 6 W is the weak mixing angle, and the functions F^ 3 , F^ 1 
and F% J are form factors that involve the Yukawa couplings and loop functions 
arising from the 7, Z penguins and box diagrams that determine the the full pro- 
cess (see appendix A for a compilation of these expressions and ref. [21] for their 
derivation). 

Following the same numerical procedure as in the U —, ^7 case and using the 
form factors given in the appendix we evaluate the fi + —, e + e~e~, r + —, e + e~e~ 
and r + — > fi + fi~fi~ decay branching ratios for both, the normal and inverted light 
neutrino mass spectra. We find that r + — > e + e~e~ and r + — > fi + fi~ fi~ processes 
are always below ~ 10 -9 (due to the constraint enforced by the seesaw condition 
when I Ai| > 10" 1 ), so in fig. 2 we only display the results for /z + — > e + e~e~ . We 
observe that the branching ratio can saturate the current experimental bound for RH 
neutrino masses M N < 0.1 TeV, 1 TeV, 10 TeV provided |Ai| > 2 x 10~ 2 , 10" 1 , 1, 
respectively, very similar to the \i — > cy case. The results for the inverted light 
neutrino mass hierarchy are again very similar and consistent with these values. As 
can be seen by comparing figs. 1 and 2, with the sensitivities of the planned future 
experiments this process has the potential to probe considerably larger values of 



3 The proposed Mu3e experiment at PSI aims for a sensitivity of 10 15 in its first phase and 
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M N [GeW] M,v[GeV] 

Figure 2: Decay branching ratio BR(fi~~ — > e~e + e~) normalized to |Ai| 4 for nor- 
mal light neutrino mass spectrum as a function of common RH neutrino mass (left 
hand side plot). The upper horizontal dashed line indicates the current bound on 
the /i~ e + e~e~ rate placed by the SINDRUM experiment [22], whereas the lower 
dotted one illustrates prospective future experimental sensitivities of the Mu3e experi- 
ment [24]- The corresponding bounds on |Ai| from the present [22] (upper gray band) 
and prospective future [24] (lower brown band) experimental searches are shown in 
the right hand side plot. The widths of the bands are due to the uncertainties in the 
neutrino mass matrix parameters. The results for the inverted neutrino spectrum 
are very similar and are thus not shown separately. 

the RH neutrino masses (compared with fi — > e^), reaching RH neutrino mass scales 
in excess of (9(10 5 GeV) for |Ai| ~ 1. Finally we note that due to the strong |Ai| 
dependence, values of |Ai| below 1CT 3 are not expected to yield observable rates at 
near future experimental facilities even for RH neutrino masses of the order 100 GeV. 



3.3 fi — e conversion in nuclei 

Competitive lepton flavor violation constraints can also be obtained from searches for 
fi—e conversion in nuclei. Currently the strongest bounds on BR^g = r convers i OI1 /T captu 

(TO 

were set by the SINDRUM collaboration from experiments on titanium with BR/t e < 
4.3 x HT 12 [25] and gold target setting BR^ u) < 7 x 10~ 13 [26], both at 90%CL. 
The — e conversion bounds are expected to be further improved in the future by 
several orders of magnitude. According to proposals [27] and [28, 29], one can expect 
a sensitivity of 10~ 16 or even 10~ 18 by the PRISM/PRIME experiment. 

To get the constraint in the fi — e channel from these experiments, one needs 
to compute the relevant transition matrix elements in different nuclei. A detailed 
numerical calculation has been carried out by [30] and we use their formula in eq. (14) 
to calculate the desired conversion rates. They receive one-loop contributions from 
photonic penguins contributing to both effective dipole (An) and vector (g^P) 

1CT 16 in its second phase [24]. 



10 



Nucleus 


D[mT\ 


y{p) 


rrifi J 


V^[ml /2 


TcapturctlO S 




Tig 


0.0864 


0.0396 


0.0468 


2.59 


Au^ 7 


0.189 


0.0974 


0.146 


13.07 



Table 1: Data taken from Tables I and VIII of [30]. 



couplings, as well as Z penguins and W box diagrams (these only contribute to 
Slv^)- Using the notation of [30] we thus have 



conversion 



2G F 



A> R D + (2g<$ + g$)vV> 



(<*)• 



(23) 



where G F is the standard model Fermi coupling constant and Ar, g^f 1 are found 
to be (Q M = 2/3, -1/3): 



.4 



V2 



aw 



(u) 

9lv 



(d) 
9lv 



8G F M W 8tt 



y/2. 



a 



w 



SG F M W 



Fr + 4F 



/i3e(l) 



y/2i 



a 



w 



8G F M^ 



Box 



r Z r Box 



- ±QuS 2 w (Ff - Ff ) 
+ AQ d s 2 w (F? - ) 



(24) 
(25) 
(26) 



Following the same numerical procedure as in sections 3.1 and 3.2 we evaluate the 
resulting jj, — e conversion branching ratios. Since both Ti and Au processes feature 
the same flavor structure, the differences between them are entirely determined by 
the numerical factors quoted in table 1. The Ti parameters entering in the conversion 
rate are on average a factor ~ 2.5 smaller than the ones for Au, whereas the capture 
rates differ by a factor ~ 5. Accordingly the difference between these branching 
ratios is a factor of ~ 2. Due to its more stringent experimental upper bound we 
thus display only the results for Au in fig. 3 for the case of the normal light neutrino 
mass spectrum (the differences with the inverted mass spectrum case are again tiny). 
For M N ~ 1 TeV the pieces proportional to V^'^ in eq. (23) cancel, so the /i — e 
conversion rate around this RH neutrino mass value is mainly controlled by the 
dipole contribution, A* R D, which is roughly two orders of magnitude smaller. The 
dip in fig. 3 is due to this cancellation. 

From fig. 3 it can be seen that the current experimental bound on this process 
imposes a constraint on the RH neutrino mass (as a function of |Ai|), which is 
roughly a factor of 2 stronger compared to bounds from \i — > and fi — > 3e (except 
for the region around Mjy ~ 1 TeV, as explained above). Furthermore, given the 
expected future sensitivities, ^"Au^g 7 — > e~Au7g (and jj'Ti^l — > e'Ti^) could 
probe RH neutrino masses up to C(10 3 TeV) , far above the values accessible in 
ji — > e7 and /i~ — > e~e + e~ , and thus constitutes the primary search channel for such 
scenarios of heavy RH neutrinos. 
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M,v[GeV] M N \GeV] 

Figure 3: Branching ratio BR([i~ Auyg 7 — > e~ hxx^ 1 ) normalized to |Ai| 4 as a function 
of the common RH neutrino mass for the light neutrino normal mass spectrum (left 
hand side plot). The upper horizontal dashed line indicates the current bound on this 
process settled by the SINDRUM experiment [26], whereas the lower one indicates 
future experimental sensitivities of the PRISM/PRIME experiment [27, 28, 29}. The 
corresponding bounds on |Ai| from the present [26] (upper gray band) and prospective 
future [27, 28, 29] (lower brown band) experimental searches are shown in the right 
hand side plot. The widths of the bands are due to the uncertainties in the neutrino 
mass matrix parameters. The results for the inverted neutrino spectrum are very 
similar and are thus not shown separately. 

4 Primordial lepton asymmetries 

We now turn to the issue of primordial lepton asymmetries and the related dynamics 
of the RH neutrinos. As we have discussed, different -R-charge assignments allow to 
define two types of models of which type B may yield sizable charged lepton flavor 
violating decays. For these effects to take place RH neutrino masses at or below 
the TeV scale as well as Yukawa couplings of order 1CT 2 or larger are needed. The 
washouts induced by such couplings and in this mass range are so large that any 
lepton asymmetry generated via the out-of-equilibrium decays of the RH neutrinos 
will always yield a baryon asymmetry much smaller than the observed one [31] 6 . 

Either producing a baryon asymmetry consistent with the observed value or not 
erasing a preexisting one via the dynamics of the RH neutrino states (in case the 
RH neutrinos are still light and the resonant condition M^ 2 — Mfj x ~ is not 
satisfied) requires small Yukawa couplings, thus rendering charged lepton flavor vio- 
lating decay branching ratios negligibly small. The phenomenological requirements 
of sizable charged lepton flavor violating effects and the generation of a B — L asym- 
metry (or of not erasing a preexisting one) are therefore mutually exclusive. Since 
these requirements cover non-overlapped regions in parameter space they are from 
that point of view complementary. 

6 In models with a slightly broken lepton number the washout is tiny, as it is determined by the 
amount of lepton number violation [13]. In our case since lepton number is broken even in the 
U(X)r symmetric phase the washouts are dominated — as usual — by A^ 12 inverse decays. 
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The generation of a B — L asymmetry in the type B models discussed here follows 
quite closely the analysis done in ref. [32]. Thus, we do not discuss this issue here 
and instead study the constraints on parameter space derived from the condition of 
not erasing an assumed preexisting B — L asymmetry. Note that in type-I seesaw 
models with flavor symmetries in the lepton sector, as for example in MLFV models, 
the CP violating asymmetry in RH neutrino decays vanishes in the limit of exact 
flavor symmetry [33]. However, since in type B models the MLFV hypothesis is a 
consequence of the intrinsic structure of the model this does not happen. 

In order to quantify these effects from now on we focus on the normal hierarchical 
light neutrino spectrum. Results for the inverted hierarchical case resemble quite 
closely the ones reported here. We start by recalling that the washouts induced 
by both RH neutrino states (at T ~ M) on any primordial B — L asymmetry are 
determined by the following set of kinetic equations: 

f'^E^W. (27) 

Here Yx = (nx ~ n x) I s (where nx is the number density of particle X and s is the 
entropy density), z = M/T and A» = B/3 — Li with Li = 2£i + e^. The function K\ 
is the modified Bessel function of the first type and the flavor coupling matrix 
is determined by the chemical equilibrium conditions imposed by the reactions that 
at the relevant temperature regime (T ~ M) are in thermal equilibrium [34]. The 
parameter /c i; that determines the strength of the flavored washouts, is given by 

fhi _ . . o 

Ki= — where m, = 2— Aa . (28) 

M 

The factor m+ ~ 1.1 x 10~ 12 GeV. Note that in the basis in which the RH Majorana 
neutrino mass matrix is diagonal couple to the lepton doublets with strength 
Aa, the factor 2 in fhi is due to this fact. 

According to the parametrization in eq. (13) the Ki parameters can be written 

as 



l + pU: 3 + y^~pU; 2 . (29) 



v 2 |Ai| 2 ,c |2 _ v 2 |Ai| s 

Ki — An — — T~r 

M m+ M 

Thus, after fixing low-energy observables the values of the parameters Ki depend only 
on M and |Ai | . Fig. 4 (left hand side plot) shows an example for the values of K e ^ (the 
k t is smaller than k^ by less than a factor 10) obtained by enforcing neutrino data 
to lie within their 2a experimental ranges [1] and fixing for concreteness |Ai| = 10~ 5 . 
As can be seen, if the preexisting asymmetry is sufficiently large even in the case of 
light RH neutrinos a sizable asymmetry in the electron flavor could be stored. 

An estimation of the N± t 2 washout effects can be easily done in the one-flavor 
approximation by taking = I in eq. (27). The resulting equation can be ana- 
lytically integrated yielding the following result for the final baryon asymmetry: 

Y AB = ^Y^_ L e~^ 8 . (30) 
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M N [GeV] M N [GeV] 



Figure 4: Left hand side plot: washout factors for muon and electron lepton flavors 
as a function of the common RH neutrino mass in the case of a normal hierarchical 
spectrum. Right hand side plot: \Y& B \ as a function of the common RH neutrino 
mass for several values of the assumed primordial B — L asymmetry. The solid 
(black) horizontal line indicates the observed value of the baryon asymmetry. See 
the text for more details. 



From this equation a parametric relation between the relevant parameters can be 
calculated, namely 

|Ai| 2 



M 




(31) 



thus fixing Fa s to its central value (Ya b = 8.75 x 10 11 [31]) and taking Y^ C 
[1(T 8 , 10~ 2 ] it turns out that as long as |Ai| 2 /M C [1, 5] x 1CT 16 GeV -1 a primordial 
asymmetry may always survive the related washouts and yield a value consistent 
with the observed one. 

A precise treatment, however, requires the inclusion of flavor. In the mass 
range we are interested in ([10 3 ,10 6 ] GeV) all the standard model Yukawa pro- 
cesses (quarks and leptons) are in thermodynamical equilibrium [34]. Neglecting 
order one spectator processes, the kinetic eqs. (27) consist of three coupled differ- 
ential equations accounting for the evolution of the A TiMj6 asymmetries. Defining 
the asymmetry vector Y^ = (Ya t , Ya, Y^J the system of coupled equations can be 
arranged in a single equation 

^--CT^^'M**- (32) 

where C\f = |Aji| 2 C^ and the matrix C^ e \ at this stage, is given by [34] 



'221 -16 

C (/) = — | -16 221 -16 ) . (33) 
-16 -16 



711 
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By rotating the asymmetry vector in the direction in which becomes diagonal 
(1a' = PVa) the system of equations can be decoupled and thus solved analytically 
for 1a' as in the unflavored regime: 

i Y *=-^M &£ * YA ' Kliz)z3 w,th p&,)p - 1 = 6 Z < 34 > 

The solution reads 

11, = V™ e~^ 8 , (35) 

where the Cj's (i = r, /i, e) are the eigenvalues of the matrix C^. The final baryon 
asymmetry in this case is therefore given by 

y a , = | E ^ = | £ e - 3 ""'- /8 • < 36 > 

j=r,fi,e j,i=r,fi,e 

In order to illustrate the effects of the Ni % related washouts on a preexisting B — L 
asymmetry we fix the light neutrino mixing angles and the atmospheric and solar 
scales to their best fit point values [1], 5 = n/2, <p = and again |Ai| = 10~ 5 . 
Assuming the same primordial Aj asymmetries in each flavor, varying them from 
10~ 8 — 10~ 2 , and using eq. (36) we calculate the resulting Y& B asymmetry The 
results are displayed in fig. 4 (right hand side plot). It can be seen that for the set 
of parameters chosen a Y"a s m the observed range can always be obtained. 



5 Conclusions 

Besides the global total lepton number U(1)l the canonical seesaw mechanism also 
breaks a global U (1)r symmetry respected by the kinetic and gauge terms in the SM 
Lagrangian. In the context of MLFV models, this U(1)r can be identified with global 
phase rotations of the charged lepton electroweak singlets e or RH neutrinos N. In 
this paper we have explored the implications of a slightly broken U(1)r symmetry 
in the context of minimal seesaw setups (with two RH neutrinos). We have shown 
that depending on the i?-charge assignments two classes of generic models can be 
identified: (type A) models where the small breaking of U(1)r allows to decouple 
the lepton number breaking scale from the RH neutrino mass scale [8]; (type B) 
models where the parameters that slightly break the U(1)r induce a suppression in 
the light neutrino mass matrix. 

We have studied the implications of these models for charged lepton flavor vio- 
lating decays. We found that in type A models the decoupling of the RH neutrino 
masses from the lepton number breaking scale implies also a suppression of the cor- 
responding Yukawa couplings, thus leading to non-observable charged lepton flavor 
violating effects. Type B models realize the MLFV hypothesis in the sense that 
due to the structure of the light neutrino mass matrix their flavor patterns are — up 
to normalization factors — entirely determined by low-energy neutrino observables. 



15 



Moreover, the suppression induced by the slightly broken U(1)r on the neutrino mass 
matrix allows large Yukawa couplings and TeV RH neutrino masses, and thus po- 
tentially large flavor violating \i processes. We have studied the /i — > ej, /i — > 3e and 
fi — e conversion in nuclei for normal and inverted neutrino mass spectra, finding that 
the three processes have branching ratios accessible in present experiments as long 
as the relevant overall Yukawa normalization factor is larger than ~ 1CT 2 , 10 -1 , 1 
and the RH neutrino masses are below ~ 0.1 TeV, 1 TeV, 10 TeV, respectively. For 
heavier RH neutrinos fi — > ej is below prospective future sensitivities while fi — > 3e 
and fi — e conversion in nuclei would remain observable, up to ~ 100 TeV and 
Mjy ~ 10 3 TeV respectively. On the other hand in both type A and B models, RH 
neutrino contributions to LFV tau lepton decays are restricted below the present 
and near future experimental sensitivities. 

Sizable fi flavor violating decays require large Yukawa couplings and light RH 
neutrinos. These values imply large RH neutrino inverse decay effects, that render 
the dynamics of these states incompatible with either the generation of a B — L 
asymmetry (consistent with the observed B asymmetry) or with the preservation 
of a preexisting one. Accordingly, sizable lepton flavor processes and small RH 
neutrino inverse decay effects are phenomenological requirements that cover non- 
overlapping regions of parameter space, from that point of view the analysis of both 
of them turns out to be complementary. In the low mass range (M < 10 6 GeV), 
instead of studying the generation of a B — L asymmetry via resonant leptogenesis, 
we have considered the influence of the RH neutrino dynamics on a primordial 
B — L asymmetry. We have demonstrated that a preexisting asymmetry yielding 
the observed B asymmetry can survive the RH neutrino related washouts provided 
the overall Yukawa coupling normalization is below ~ 10~ 5 . 
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A Formulas for k — > Ijj and /~ — > Ijlflj processes 

In this appendix we summarize the formulas we use for the calculation of the charged 
lepton flavor violating decays discussed in sections 3.1 and 3.2. The results presented 
here were extracted from ref. [21] and adapted to our notation. In what follows the 
parameters r a 's are defined according to r a = M^/M^ . 

The process /r — > Ijl^lj is determined by 7, and Z penguins and box diagrams 
(for the full set of Feynman diagrams see ref. [21]). The 7 penguin contribution can 
be split in two pieces corresponding to the photon being either on-shell or off-shell. 
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For the on-shell piece, the one that determines the U — > Ij^y process, we have 

2 



^ = -(A-G 7 .A%, 

G 7( r a) = jt~i °„ \a ( 2 + 3r ' a ~~ 6r a + r l + 6r a logr a ) , 



4(1 -r a y 
whereas for the off-shell photon piece 



(37) 
(38) 

(39) 



[7 - 8r a -Ur 2 a + \2r\ - (2 - 20r a + 2Ar 2 a ) \ogr a ] . (40) 



12(1 -r a y 

The Z penguin contribution can be split in two parts, namely 



jpklj JphhW _|_ JpL 



i(2) 



where the first piece can be written as 

ihW _ 2 



F z (r a ) 



(l-r a + logr a ) 



2(l-r a ) : 



1-r 



log r a , 



while the second contribution according to 



plilj(2) 



9' 



V(A) 
Z 

Gf(r a ,r b ) 
Gf{r a ,r b ) 
Gf{r a ,r b ) 



Gf\r a ,r b ) = {X ■X) ab G ( z A \r a ,r b ) with A = 2,3,4, 

logr -- — logr 6 



2(r a -r b ) \l-r a 
r a r b 

logr a , 



1 - n 



2(1 -r a ) 



log r b , 



2(1 -n) 

Hz(r a ,r b ) = (\ T -\*) ab H z (r a ,r b ), 

r b (l - 4r a ) 



H z (r a ,r b ) = 



4(r a - r b ) 



l — r„ 



\ogr a 



r a (l - 4r b 
1 - r 6 



log r b 



(41) 

(42) 

(43) 
(44) 

(45) 
(46) 
(47) 

(48) 
(49) 
(50) 
(51) 



Note that due to the constraint implied by the SU(3)e+N flavor symmetry the off- 
diagonal elements of the matrices G^\r a , r b ) and H^\r a ,r b ) vanish. 
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The box diagram contributions can be split in three parts as follows 

£ 3 



^Box 



^li3lj (A) 
Box ' 



(52) 



,4=1,2,3 



For the first part we have 



-^Box 



^Box( r a) 



2r, 



(1 - r a ) 



^— - (1 - r a + r a \ogr a ) 



For the second is given by 
4 



(2) 
^Box 



F^(r a ,r b )(j)^X* a F^2(r a ,r b )X jb with A = 2, 3 



^bI(^) 



r a r 6 



4(r a - r 6 ) 



1 - 4r a (2 - r b ) 1 - 4r b (2 - r a ) 

log r a log r b 



(l-r a )(l-r 6 ) 
?'6 



(l-r a ) 2 
(7 - 4r a r fe ) 

-(1 - r b + logr fe ) + 



(1 - nf 



(53) 
(54) 

(55) 
(56) 

(57) 



;i-r a ) 2 



1 - r a + r a \ogr a ) 



(58) 



where in F^^(r a ,rb)(J) 110 summation over the indices a, b is performed. Finally, 
the third term in (52) can be written as 



T-i'i3ij(3) / -\ 
^Box UJ 



A ■ Gbo X (j) ■ At 



G , Box(r a ,r fe )(j) = \ ja G B ox( r aj r h) Xj b 



(59) 
(60) 



G B ax(r a ,r b ) 



r a - n 



r a [1 - 2r 6 (l - 2r )] . r 6 [1 - 2r«(l - 2r b )] . 

- log r a i log r fe 



(l-r a ) 2 
;i + 2r a r b ) 



(61) 



(1 - r a )(l - r fe 

where, again, in G^(r a , r&)(j) no summation over the indices a and b is performed 
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